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Abstract. We report on a simultaneous measurement of the top quark, W boson, and 
neutrino masses in ti dilepton decays with 4.98fb _1 of data from the CMS experiment. The 
analysis is based on endpoint determinations in kinematic distributions. In addition to the 
unconstrained fit for three masses, the neutrino and W boson masses can be constrained 
to standard values; in the maximally constrained fit, the top quark mass is found to be 
M t = 173.9 ± 0.9 (stat) +\l (syst) GeV. 



1. Introduction 

In many proposed extensions of the Standard Model, there is a new Z2 symmetry, like R parity 
in Supersymmetry (SUSY), under which the new physics is odd and the known particles are 
even. New physics particles are thus produced in pairs and decays must terminate in a stable 
particle. Such particles are candidates for dark matter and escape direct observation in the 
CMS detector [lj. We have developed a model independent method to simultaneously extract 
all unknown masses in a decay chain. With the production of pairs of neutrinos, dileptonic ti 
resembles the generic new physics topology. Thus, we test this method in data by simultaneously 
computing the top quark, W boson and neutrino masses. Going beyond the general study of 
the technique, we also fix the W and neutrino masses to their known values and perform a 
measurement of the top quark mass. This measurement is competitive in the dilepton channel 
and does not rely on Monte Carlo simulation. More details are available in Ref. [2]. 

2. Kinematic Variables - Mt2 

To combat the problem of underconstrained kinematics, we look at extreme kinematic features 
of distributions. In particular, endpoints of distributions are determined and related to the 
underlying masses. This analysis uses the extension of transverse mass called Mt2 |3]- M?2 
is the minimum parent mass consistent with observed kinematics. We have three unknowns, 
M„, Mw, and M t , so we partition the ti decay in order to create three mass variables |4j. Two 
Mt2 subsystem variables are combined with the bl invariant mass to simultaneously extract the 
three masses. Figure [T] details the partitioning of the event topology for the Mt2 variables. In 
general, the resulting variables depend on the magnitude of the transverse upstream momentum, 
Pt [5]. To remove this dependence, we project all momenta perpendicular to Pt and compute 
Mt2±- From [1], we have analytic formulae that relate distribution endpoints and underlying 
masses. 
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Figure 1. We partition the ti decay into subsystems. For a third variable, we use the bl 
invariant mass, which is in many cases the same as the third possible Mt2 variable we could 
form. 



3. 2011 Data (4.98/fb) 

Monte Carlo (MC) is overlaid on the 2011 data to illustrate the approximate signal and 
background content of the distributions in Figure [2j There is good agreement, but the signal 
MC is not used in the reported measurement. From the known values of Mt, Myy, and M u the 
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expected endpoints of Mj, 2 ± , Mj, 2 ± , and mu are 80.4, 173, and 152.6 GeV (respectively). 
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Figure 2. Monte Carlo overlaid on the 2011 data for the three measured kinematic distributions 
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j I j an d Mbi- The MC is luminosity normalized to the data with NLO cross sections. 



4. Fitting Strategy 

We perform an unbinned maximum likelihood fit that takes into account event-by-event 
information such as resolution function shapes. For M = (M^, Mw, Mtj^ and data vector u 
the likelihood is given by £(M) = UaU^=i ents ^i( u i\ M ) for a G i M T2 > M fh M bi} and 



C? (xi\x max ) =a j S a (y\x^ x )TZ^(x l - y)dy + (1 - a)B a ( Xi ), (1) 

where a is the signal normalization, 1Z is a function that incorporates event-by-event resolutions, 
B is the background estimation and S is the signal shape. This last quantity is approximated 
as linear decent into the endpoint in a region near the a priori kinematic maximum. 



1 We are only sensitive to the square of the neutrino mass. 



5. Background Estimation 

The dominant source of background is combinatorial. In particular, the largest contribution to 
this background is from b jet selection, which is not 100% pure. Non ti events do not obey the 
same bounds as ti. Figure [3] shows the composition of events in the invariant mass distribution 
and makes clear that the mis-tag background dominates the events beyond the endpoint (153 
GeV) . We model this dominant background by using a control sample in which we have inverted 
our b tagging. We then use a non-parametric estimation (adaptive kernel density estimation) of 
the shape which enters as B a in Eq. 6. 
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Figure 3. Monte Carlo composition of the My. Note that the dominant source of events beyond 
the endpoint are due to signal events in which the wrong jets were chosen as the b jets. 



6. Uncertainty Estimation 

The usual methods for statistical uncertainty estimation are not valid because they do not 
account for fluctuations in the background shape. We thus use bootstrapping to estimate our 
statistical uncertainty. By recomputing the fit to the control region for each bootstrap sample, 
we can take into account background fluctuations by taking the central value and width of the 
distribution of masses evaluated from the ensemble of bootstrapped samples. 

We estimate our systematic uncertainties by inserting variations into our fit and re-evaluating 
the result. Table [I] summarizes the contributions to the total uncertainty. Except for color 
reconnection, all the variations are applied to the data. 



Systematic 


Uncertainty (GeV) 


Systematic 


Uncertainty (GeV) 


Jet Energy Scale 


+0.5 -1.4 


Fit Range 


±0.6 


Jet Energy Resolution 


±0.5 


Background Modeling 


±0.5 


Pileup 


neligible 


Color Reconnection 


±0.6 



Table 1. Systematic Uncertainties. Except for color reconnection, the variations used in 
determining the uncertainties are applied to data. The total uncertainty is +1.2, —1.8 GeV. 



7. Validation 

As a cross check of our analysis, we perform our fit on many pseudo-experiments (PE) of 
simulated data. For each PE, define Pull = (mt rae — m meas )/a meas . If our procedure is unbiased 
and correctly produces the uncertainty, then the distribution of the pull should be a standard 



normal. We can also probe potential bias by repeating this exercise for various values of mt rue . 
As shown in Figure [4j both of these tests are consistent with an unbiased measurement with 
correct statistical uncertainty. 
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Figure 4. The left plot shows the pull distribution, with a mean of 0.15 ± 0.19 and a standard 
deviation of 0.92 ± 0.06. The right plot shows the spread in measured top masses as a function 
of input top mass. The best fit line is shown in red (slope 0.91 ± 0.4, intercept 0.34 ± 0.3) and 
the line of unit slope is in blue. The blue line agrees with the fit at x 2 / n df = 10.7/9. 



8. Results 

We present three sets of results, corresponding to three levels of constraints on the masses 
of the neutrino and W boson. The set of masses are shown in Table |U In the case where 
both the neutrino and W boson masses are constrained, we compute our most precise top 
mass measurement: Mt = 173.9 db 0.9 (stat) jVg (syst) GeV. Due to its minimal dependence 
on MC, this result has many complementary (systematic) uncertainties to other measurements 
and so can help reduce the world uncertainty. However, in practice this is non-trivial as there 
is ambiguity on the order of Aqcd as to what this analysis measures and what traditional MC 
calibrated approaches compute. As the precision of top mass measurements improves, the top 
quark community will need to think carefully about the interpretation of the quantities being 
measured. 





Constraint 


Fit Quantity 


None 


m u = 




m w = 80.4 


ml (GeV 2 ) 


-556 ± 473 ± 600 





(0) 


m\y (GeV) 


72 ± 7 ± 9 


80.7 ± 1.1 ± 1 


(80.4) 


mt (GeV) 


163 ± 10 ± 11 


174.0 ±0.9 ±2 


173.9 ± 0.9±J;| 



Table 2. Results for the various levels of constraints. Uncertainties are first statistical and then 
systematic. 
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Figure 5. The fit plot showing the data and fit function for the case in which the neutrino 
and W boson masses are fixed. The displayed fit and data are from a particular bootstrapped 
sample. The blue line is the background estimation from the control sample, the green line is 
the signal fit function and the red is their sum. The top right of the mu plot shows the endpoint 
region to illustrate the fit to the data. 
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